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Abstract 

The concept of localization in Fock space is extended to the study of the 
many particle excitation statistics of interacting electrons in a two dimen- 
sional quantum dot. In addition, a finite size scaling hypothesis for Fock 
space localization, in which the excitation energy replaces the system size, 
is developed and tested by analyzing the spectral properties of the quantum 
dot. This scaling hypothesis, modeled after the usual Anderson transition 
scaling, fits the numerical data obtained for the interacting states in the dot. 
It therefore attests to the relevance of the Fock space localization scenario to 

the description of many particle excitation properties. 
PACS numbers: 72.15.Rn,71.55.Jv,73.20.Dx 
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The concept of localization in Fock space has been recently utilized [0] to explain the 
transition in the width of excited states measured in tunneling conductance experiments for 
quantum dots [0] . Due to electron-electron interactions the quasi-particle states are coupled, 
which leads to their finite width. At a critical value of interaction strength and excitation 
energy the width becomes essentially infinite and the quasi-particle states can no longer be 
resolved. 

In this letter we show that Localization in Fock space explains also the recently reported 
10-^ transition in the level statistics of excited many particle states as a function of the 
interaction strength. For small values of the interaction the many-particle states are localized 
in the non-interacting Fock space, namely each interacting state is composed out of a small 
number of non-interacting eigenvectors. Above some critical value they are extended, i.e., 
composed out of a large number of non-interacting eigenstates. This leads to a transition 
in the statistical properties of the energies from Poisson to Wigner statistics, which as in 
Anderson localization is a true second order phase transition. We will demonstrate this by 
constructing a finite-energy scaling theory which plays the role of finite size scaling familiar 
in standard localization theory. 

A typical gray-scale map of the statistical properties in the space of interaction strength 
and excitation energy is presented in Fig. [l|. at low excitation energies the statistics re- 
mains Wigner over a large range of interaction strength. This will be discussed in detail 
elsewhere. A transition from Poisson to Wigner and back to Poisson as the interaction 
strength is increased is evident for intermediate excitation energies. The second (Wigner- 
Poisson) transition is the result of Wigner crystallization . In this letter we will show that 
the first transition has the characteristics of a localization transition in Fock space. 

A localization transition is characterized by a single parameter scaling behavior |^, which 
results in only three possible types of statistics in the thermodynamic limit: Poisson in the 
localized regime, Wigner in the extended one and a possible third type of statistics at the 
transition point As shown by Shklovskii et. al. ||^, for a finite system the transition 

between the two types of statistics is gradual and may be characterized by a finite size 
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scaling behavior, which can be used to locate the critical point and the critical exponent. 

For the Anderson localization transition the scaling parameter is the ratio of the disor- 
der dependent localization length to the linear dimension of the system. In the Fock space 
localization transition we show that the relevant scaling variable is the ratio of the local- 
ization length (which depends on the interaction strength and the two-particle density) to 
the linear "dimension" of the system in Fock space, which depends mainly on the excitation 
energy. The latter is surprising, since one would expect the linear dimension of the system 
to depend on the number of particles and the number of sites. As has been demonstrated 
in Ref. P], no clear finite size scaling behavior is observed as function of the number of sites 
or particles. This is the result of the unusual "geometry" of the interaction coupling which 
leads to a Cayley tree structure in Fock space where the number of generations depends on 
the excitation energy 

Let us begin by defining the Fock space. The many-electron second quantization Hamil- 
tonian: 

H = Ho + Hint = + H Ui^j^k,lc]c\ckCi, (l) 

where c| is the creation operator of particle in the i-th single electron state, is its energy, 
and Uij^k,i is the interaction matrix element, has a set of eigenvectors I^E'^v) ^^nd eigenvalues 
E-j^, where N is the number of electrons and j is an index ordering the states by ascending 
eigenvalues. There are M = (ff) many-electron states, where S is the number of single 
electron states. 

For the non-interacting {Hint = 0) case the ground-state eigenvector and eigenvalue are 

|$°.)=4...cl|0) 

N 

E% = Y.^^i (2) 
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here |0) is the vacuum state. An excitation composed of m electron- holes pairs may be 
represented as: 

\^N I — (-a2m ■ ■ ■ ^am+&am ■ ■ ■ (-ax\^ N I 



m 2m 
1=1 i=m+l 

These excitations may also be numbered in ascending order |$7v) according to their energy. 

For the interacting {Hint 7^ 0) case one may write the eigenvectors 1"^^) ^^^^ their 
excitation energy Ej^ — (again arranged in ascending order) in the following way: 

N M 
m=l ai...a2m k=l 



e 



Eh - K- (4) 



An interesting question is how are these interacting eigenstates composed out of the non- 
interacting Fock space eigenvectors. It is natural to expect that for small values of interaction 
and excitation energy an interacting eigenstate will t)^ composed of a small number 

of non- interacting eigenvectors \^%) where k is in the vicinity of j, while for large values of 
the interaction |\E']v) will be composed of many non-interacting eigenvectors. Altshuler et. 
al. have termed this transition "localization in Fock space" and predicted that it should 
exhibit the characteristics of an Anderson transition on a Cayley tree. 

We will use two measures for the Anderson transition: (i) the energy level statistics, 
and (ii) the inverse participation ratio in Fock space. A convenient way to characterize 
the change in the statistics of a system proposed by Shklovskii et. al. is to study the 
parameter 7 defined as 

e — e 

where P{s) is the distribution of the normalized level spacings s = E^ — Ej^^ / {Elf — Ej^'^), 
where (. . .) denotes an average over different realizations of disorder. For an infinite system 
7 changes sharply from 7 = 1 in the localized regime to 7 = in the extended regime. For 
a finite system the change is gradual. One can then use a finite size scaling hypothesis to 
identify the transition point and its critical indices [^]. The inverse participation ratio P, 
defined as 

M 

p = T.\M\ (6) 

k=l 



is expected to change in the Anderson transition for an infinite system from P which is a 
function of the locahzation length in the locahzed regime to P = in the extended regime. 
The argument for a locahzation transition as one increases the interaction strength follows 



the usual Anderson transition argument |T^. The strength of the hopping term V coupling 
neighboring sites is compared with the inverse of the density of states of the neighboring 
site, Ku, where u = l/W, W being the strength of disorder, and K is the connectivity. A 
localization transition will occur at a critical value 

KV 

Zc = KuV = ^^ 1. (7) 

As noted by Altshuler et. al. [|^, a similar argument may be applied to Fock space localiza- 
tion. The two-body interaction couples non-interacting eigenvectors which are different by 
up to 2 electron-hole pairs. Assuming that the main coupling is between states differing by 
a single electron-hole pair [||], the average strength of coupling 

U = (($^--"-|if,„,|<|.^-"-))„,,.,^, (8) 

where (. . denotes averaging over all possible electron- hole pairs. The density of such 

coupled states at excitation energy e is denoted by Kh'2{e). Thus 

~ Ku2{e)U. (9) 

This argument is similar to the one used by Imry in the context of two particle state de- 
localization [Tl|J^. A more careful consideration of the geometry of the connected states 



reveals that it is similar to a Cayley tree for which K should be replaced by In if [jT2 
Since the effective connectivity K g where g is the dimensionless conductance, this is 
a correction of order of \ag. Altshuler et. al. [|1| concluded that due to the Cayley tree 
structure, an intermediate region Z^-jK > i'2{e)U > Z^/KhiK between the localized region 
i'2{e)U < Zc/ K In K and the extended region Zc/K < h'2{s)U should exist, which is non- 
ergodic. Using a non-linear sigma model for the Cayley tree Mirlin and Fyodorov found no 



intermediate region |l3 



Here we would like to extend this Fock space localization picture to include one of the 
most useful concepts in the Anderson localization picture, namely that of finite size scaling. 
In the usual Anderson picture 'y{KV/W, L) = f[L/C,{KV/W)] P where L is the sample 
size, ^{KV/W) is the localization length and / is a scaling function. Near the critical value 
it is expected to behave as -f{KV/W, L) = 7(Z„ L) + C[Z, - {KV/W)]L^/^, where C is a 
constant and 6 is the critical exponent (usually denoted by u). From the above discussion it 
is natural to assume that Ki'2{e)U will replace KV/W for the Fock space localization. The 
size of the system, i.e., the number of generations in the Cayley tree, is equivalent to the 
number of electron-hole pairs which can be generated at a given excitation energy Thus 



L ~ Umax ~ V^/^) where A is the single electron level spacing. Hence we expect 

^{Kiy2{e)U,e) ~ ^{Z,,e) + C [z, - Ku2{e)u) e'/^'. (10) 

For P{Ki'2{e)U,e), the situation is different. According to the Breit-Wigner formula 
14,15| in the localized regime the inverse participation ratio depends on U and Kh'2{e) but 



not on the linear dimension e, as long as C,{Ki'2{e)U) < L. Once ^{Kv2{e)U) > L, P will 
also depend on e. Since ^ oc (Z — Zc)^^ we expect that P for different values of e to coalesce 
for Z <^ Zc and to begin to fan out for Z ~ Z^. 

In order to check this finite size scaling hypothesis we have numerically calculated 7 and 
P for a specific many particle tight-binding Hamiltonian: 

H = Y. '^kjaljakj - 5I(4,i+i«fc,j + 4+i,j«fc,i) + h.c + Hint, (H) 

where u^j is the energy of a site {k,j), chosen randomly between —W/2 and W/2 with uni- 
form probability, and is a constant hopping matrix element. The interaction Hamiltonian 
is given by: 

k,j>l,p ^i,p\/'^ 

where U = e^/b and b is the lattice unit. 

We consider a 4 x 3 dot with 5 = 12 sites and = 4 electrons. The M x M Hamiltonian 
matrix is numerically diagonalized and all the eigenvectors I^I/jv) eigenvalues E^^ are 



obtained. The strength U of the interaction is varied between — 30V. The disorder 
strength is chosen to be = 5V in order to assure metalhc behavior for the non-interacting 
case. This is important in order to obtain a reasonable connectivity K g. For each value 
of f/, the results are averaged over 1000 different realizations of disorder. 7(?7, e) is directly 
calculated from the level spacing distribution for different values of U and e, while P{U,e) 
is calculated from the overlap of |^E']v) with the non-interacting eigenvectors. 

In Fig. H the average value of for an excited state 1'^'^^) (where i is chosen so 

that Ei ~ 8V in the non-interacting case) as a function of the interaction strength is plotted. 
As expected for small values of interaction |y4fe_jp is strongly peaked around k = i. Thus, 
the interacting state is composed of a small number of non-interacting states of similar 



energies, which has on the average a Breit- Wigner form [ll4| , p!5| . As the interaction increases 
the interacting state is composed of an increasing number of non-interacting states although 
there is no obvious sharp transition. Similar features appear in for a single realization, 

although the curves are noisier. 

In order to locate the transition point we use finite size scaling. First we calculate 
Ki'2{e) ~ J{e) by evaluating the number J{e) of non- interacting states in a given energy 
region which are coupled to each other by a single electron-hole transition. The numerically 
calculate J{e) is plotted in the inset in Fig. |^. For small values of e, J increases linearly, 
while near the middle of the excitation band it flattens out. Since the interaction strength 
U (X U one may replace the scaling parameter Kv2{e)U by J{e)U . 

The parameter 'y{J{e)U,e) is shown in Fig. ^ for values of 5V < e < lOV. This range 
was chosen in order to avoid the special region of low excitations and to avoid the symmetric 
region above the middle of the band (see Fig. |l]). Indications of finite size scaling behavior 
are observed. For large values of J{e)U, 7 becomes smaller as the size e increases. On 
the other hand for small values of J{e)U, 7 tends to be larger as the size increases. The 
point at which the lines intersect is the transition point. The quantity '~f{J{e)U,6) is then 
considered as a certain scaling function f{Q of the scaling variable ( = e^^'^^{Zc — J{e)U). 
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Practically, it is useful to shift the variable ( to y{() = where a and b are respectively 

the minimum and maximum values assumed by C. Evidently, y{() ranges between —1 and 
1. Then one expands f{() in a series of Chebyshev polynomials Tn[y{C)] (n = 0, 1, 2, ...Q). 
Minimization of the set of differences |/(C) ~l{J{^)U,e) \ results in the unknowns Zc, S and 
the expansion coefficients (namely, the scaling function itself). In all cases, it is sufficient 
to cut off the number of polynomials at Q = 10. The best fit is obtained for Zc = 1.3 ± 0.4 
and 5 = 1 ± 0.3. The function f{() and the original data points are plotted in the inset of 
Fig. |. 

The results for the inverse participation ratio P are presented in Fig. ^. For J{e)U < Z,, 
all lines coalesce, since no dependence on the linear dimension e is expected. For J{e)U ~ Zc 
the dependence on e is apparent and the lines corresponding to the different values of e fan 
out. The region in which the lines fan out correspond to Zc ~ 1.5, which is close to the 
values obtained for 7. It is interesting to note that 7(^c) ~ 0.6 is considerably higher than 
for the Anderson transition P]. 

In this analysis we have not clearly identified a crossover region between the localized 



and extended regions predicted in Ref. M. This may be explained by the results of Ref. [13 



although we do not see the jump in P predicted there. We speculate that the additional 
couplings between generations which are not taken into account in the pure Cayley tree 
picture tend to restore ergodicity, which is a key ingredient in the scaling picture. Only 
studies on much larger systems will answer this question. 

In summary, Fock space localization has been shown to describe the features of the exci- 
tation statistics of many particle systems. We have proposed a finite size scaling hypothesis 
in which the role of system size in the usual Anderson scaling picture is played here by 
the excitation energy. The scaling hypothesis was tested for a specific tight-binding many 
particle Hamiltonian and resulted in satisfactory agreement. 

We are grateful to A. Kamenev, A. Miiller-Groeling and D. Shepelyansky for useful 
discussions and correspondence. 
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FIG. 1. A gray-scale map of the parameter 7 defined in Eq. (||) for the many particle Hamilto- 
nian given in Eq. ([Tl|). Low values of 7 correspond to Wigner statistics, while values close to one 
correspond to Poisson statistics 
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FIG. 2. The non-interacting Fock space composition of an interacting eigenvalue of excitation 
energy e ~ 8V. The composition was averaged over 100 reahzations of disorder. Inset: 

the number of states coupled by a single pair transition J within an energy bin (dark histogram) 
compared with the total number of states (white histogram). 
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FIG. 4. The inverse participation ratio P or different values of interaction U and excitation 
energy e. 
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